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Abstract 

We study the applicability of Pade Approximants (PA) to estimate a " sum" of 
asymptotic series of the type appearing in QCD. We indicate that one should not 
expect PA to converge for positive values of the coupling constant and propose to 
use PA for the Borel transform of the series. If the latter has poles on the positive 
semiaxis, the Borel integral does not exist, but we point out that the Cauchy 
pricipal value integral can exist and that it represents one of the possible "sums" 
of the original series, the one that is real on the positive semiaxis. We mention 
how this method works for Bjorken sum rule, and study in detail its application 
to series appearing for the running coupling constant for the Richardson static 
QCD potential. We also indicate that the same method should work if the Borel 
transform has branchpoints on the positive semiaxis and support this claim by a 
simple numerical experiment. 

PACS numbers: 12.38 Bx, 12.38 Cy; 
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1 Introduction 



Perturbative expansions appearing in QCD and many of its models simmulating 
various interesting features of this theory have, usually, coefficients growing so fast 
that it is diffcult to obtain an increased precision by naive adding of subsequent 
terms [T] . | . 

Pade Approximants (PA) are an effective tool to "sum" an asymptotic series 
0,0], however they may not necesserily sum it to the function that has the required 
analytic properties. We must remember that there is an infinite number of analytic 
functions having the same asymptotic expansion, and the additional requirements 
stemming from the physical meaning of the function can help us to diminish this 
abundance. 

We shall first illustrate this problem on the very simple example and show that 
the use of the Borel integral can be of help here. As we expect the perturbation 
series in QCD to have coeffcients growing like n\K n rC f2| let us consider the simplest 
series of this type - the Euler series: 

oo 

e(z) = J2nh n (1) 

n=0 

It is well known 0,0 that PA to this series converge to the function: 

E(z) = / du (2) 

Jo 1 — zu 

in the whole complex plane with the positive semiaxis excluded. The function has 
the branch point at z = and PA reconstruct the branch of this function which 
is real on the negative semiaxis and has cut along the positive one. In QCD the 
expansion parameter is usually the coupling constant and we want to have our 
function for the positive values of coupling constant - moreover we want, for real 
expansion coefficients, the function to be real, at least close to zero. Our aim is 
therefore to find a function which has not only the expansion ([!]) but is also real 
on the positive semiaxis. We want to stress that such function is again not unique, 
but we want to show how one of such functions can be constructed from ([[]). 

We first consider when the Borel transform of the series has only poles 

(Section 2) and illustrate our point in Section 3 by two examples of QCD related 
series. Then, in Section 4, we apply our considerations to the case when the Borel 
transform has branchpoints and claim that the same approach (taking the Principal 
Value of the Borel integral over Pade Approximants to the Borel transform of the 
series) would work here - we support this claim by a simple numerical example. 
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2 Poles in the Borel transform 



We consider the series: 

oo 

f(z) = £ f n z n (3) 

n=0 

where we expect that asymptotically (for n — > oo) /„ behaves like n\K n rO ' . The 
Borel transform of this series is then: 

OO OO £ 

9(z) = Y,9nZ n = Y,^ n (4) 

n=0 n=0 n - 

If all / n 's are positive (which is often the case in QCD), then we can expect that 
g(z) has a singularity at z = K (where K is the radius of convergence of (fD). In 
this case the Borel integral: 

/•OO 

= / e - u g(uz)du (5) 



is complex for z > 0, if it exists. If g(z)(z — K) is finite for z K, then this 
singularity does not spoil the convergence of ([5]), and similiar conditions can be 
imposed on the behaviour of g(z) near other singularities. 

Let us first assume that g(z) has a simple pole at z = K. In this case: 

lirn ig(z)(z -K) = Res g(z) < oo (6) 

z^K Z= K 

and the integral @ can be written , for xe(0, oo) as: 

$(x) = PV / e~ u g(ux)du ± —e~ K / x Res g(x) (7) 
Jo a; x=_ft: 

- choice of the sign depends on whether in (|]) we integrate just below or just above 
the positive semiaxis. We see that the imaginary part of $(2) has all coefficients 
of the power expansion at z — 0+ vanishing and therefore: 

S{x) = Re${x) = PV / e- u g{ux)du (8) 

Jo 

is the function we need for x > 0. 

If we are frustrated by the fact that S(x) is defined only for x > 0, and is not 
an analytic function in the complex plane, we can define it also there: 

S(z) = / e- u g(uz)du + —e~ Klz Res g(x) (9) 

JO+ie Z x=K 

Obviously we can proceed in an analogous way if g(z) has more poles on the 
positive semiaxis. Let us, for the moment, assume it has only poles there. 
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We now ask how can we construct a sequence of approximations to £(z) from 
(^|). We have already pointed out that the direct application of PA to (0) will most 
probably fail, because we can expect that PA will diverge on the positive semiaxis. 
However, as we know that PA are particularity efficient in approximating functions 
in domains of their meromorphy ||, we can use PA to approximate g(z) from (Q), 
i.e. we propose to construct [M/N] g (z) and approximate £(z) by: 

/•CO 

S M)N (z)=PV e- u [M/N] g (uz)du (10) 

<J 

The convergence of PA in measure in the neighbourhood of the positive semiaxis 
(where we assumed g(z) to be meromorphic, and therefore we expect this type of 
convergence) guarantees us that (|10D will converge to @. We want to point out 
that the idea of combining PA with the Borel method was proposed quite a time 
ago but in a context where the Borel transform was regular on the positive 

semiaxis. Essentially the same concept was used in ||, where the Partial Pade 
Approximant [|9] was used to predict a sum of the series for the Gell-Mann-Low 
function of the effective charge. 

Before we consider a more general situation when g(z) has also branch points 
on the positive semiaxis, let us see how the method works for two examples from 
QCD. 



3 Two QCD series with only poles on the positive 
semiaxis 



First, we consider the Bjorken sum rule. It is known [TIJ that in the large-/?o 
approximation it has only four simple poles in the Borel plane and therefore our 
prescription would be exact for all [M/N] with M < 3 and N < 4. For realistic 
values of /3 , the perturbation series has been calculated up to 3rd order and there 
exists an estimate for the next coefficient 



111. It has been demonstrated in 11 



[2/1] PA to the Borel transform of the truncated series for f(x) (we take iV 



/ 



f(x) = l-x- 3.58a; 2 - 20.22x 3 - 130x 4 



where f(x)is defined in: 



[gl(x, Q 2 ) - g?(x, Q 2 )}dx = - \ g A \ f(a s / vr) 

f! "0 



that 
= 3): 

(11) 



(12) 



exhibits a pole at y = 1.05 (y = ), while we expect poles at ±1,±2 in the 
large-/?o limit, as mentioned above. Therefore we believe that 



/(*) 



1 - zPV j™ e- u [2/l\ s {u^jz)du 



(13) 
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is a reasonable approximation to f(z) (for the exact definition of S see ||11||). We 
shall not discuss this example in more depth, but only remark that the value of 



a 



: (Q 2 = 3GeV 2 ) found this way is .371; 



.075 
.068- 



The next example, which we want to discuss in detail, is the asymptotic series 
for the running coupling constant obtained from the static QCD potential for which 
we take the Richardson potential |fL2| . We follow notations and conventions of ||13|| , 
and for more clarity we shall cite some formulae from this paper. 

The Richardson potential in momentum space is: 



V R (q) 



16ll 2 C F 



(3 q 2 ln(l + q 2 /A 2 R ) 



(14) 



and we use rif = 3, therefore /3 = 11 — 2/3rif = 9, also A R = .4 and Cp = 4/3. 
This example is very instructive because one can find the exact formula for the 
perturbative coupling in the position space: 



Po 



du 



Ajjrn 



'l U ln 2 (u 2 — 1) + TT 2 

At the same time the formula for the coupling in momentum space is: 

1 



(15) 



4rr 

Po 



ln(l + q 2 /A R ) g 2 
From this, one can calculate the (3 function: 



(16) 



Pr{<*r) = q 



.da R {q) 
dq 2 



•Itt 



- a R 



ftp 2 | A> 

'4:ir aR 4tt l + q 2 /A 2 R 



(17) 



On the other hand, one has (for any potential) the expansion: 

d 

n=0 



dai 



a v (q = l/r') 



(18) 



where r' = re lE , and 7^ is Euler's constant. / n 's are coefficients of the expansion 
of: 

(2u) 2n+1 



Itan(nu) 



-exp 



E 



-C(2n + 1) 



/I fn 



u 



(19) 



Lti 2n + l j n=0 

For details see [T^ . 

For later considerations it will be useful to observe that f(u) can also be written 
in a different form: 

1 00 1 I u 



k=l 



2k+l 
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From this formula we can immediately see that f(u) has zeros at negative integer 
values of u and simple poles at half-integer positive values. 

It is important to observe from ( |TT| ) that /3r(or) is not analytic at atR(q) = 0. 
The reason is that — > for q — > oo and therefore: 

Po 2 — 47r 

{3r(o£r) = ——a R + terms like e (21) 

4:71 

Summing up, we see that we can get for a^(l/r) the perturbation expansion when 
we neglect "higher twists": 

^m = t^){^^Jn\fn (22) 

with q — 1/r'. The direct calculation of PAs to this series confirms our expectations: 
the majority of zeros and poles of diagonal PAs [M/M] (the same is the case for 
paradiagonal sequences of the type [M ± 1/M]) lie on the positive semiaxis and 
interlace: 
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zeros 



[2/2] 
[3/3] 

[4/4] 
[5/5] 



poles 



[6/6] 



[7/7] 



1 ^ 9491 n 




1 Q1^8Q 






n 17087 

U.l 1 VJO I 




fl 1 57D9 






n rn&^9 

U.UOOOZi 




U.UOc/lO 


-1-0 


9^Q90i 




-D 99^4^i 


U.UOc/lO 


yj 


9^Q90i 


PI P1Q7&ZL 




.UjUOO 






n i Am q 


994Q4i 


fl 14^Q7 

-U. l^Oc/ 1 


n 


91 ^87i 


H 1 401 Q 




14^Q7 
-u. i^tjy / 


n 

-u 


91 ^87i 


n i ^71 1 




fl 1487Q 






U.UJ i i o 




.\JO i i u 






_n 9fiQ41 


f| 9(141 4i 


-0 944Q9 


n 
u 


1 84fi8i 


-0 9fiQ41 


-fl 9D41 4i 


_n 944Q9 


-0 
u 


1 84fi8i 

10t:U01 


VJ . ZjVJ i OO 




VJ. IOjOU 










fl 0770Q 






VJ .UlUOaJO 




fl D4fl^QQ^ 

VJ .UlUOaJO 






VJ . O iiVJ It: 


VJ. IOaiUaiI 


_n 97890 


n 


1 71 SQi 


-D S9R14 


_n 1R989i 

VJ. IOaiOaiI 


_n 97890 




1 71 3Qi 

X 1 XO.J1 


n oqoo/1 










0.09094 




0.09047 






0.05022 




0.05021 






0.03037432 




0.03037431 






-0.32989 


0.18119i 


-0.28017 





17065i 


-0.32989 


-0.18119i 


-0.28017 


-0 


17065i 


0.24307 




0.21733 






0.09256 




0.09205 






0.05154 




0.05153 






0.03165505 




0.03165504 






0.01674128 




0.01674128 







We see that for large q (a^ small), there will be difficulties in finding stable 
values for PAs, because of condensation of their zeros and poles. 

When we now go to the Borel plane we get the series which is just: 



a R {q)f{ 



ffoag(g)- 

47T 



(23) 



However we know already (see Q2DD) that f(u) has only poles on the positive semi- 
axis, and therefore we can try to sum (B2) using the formula: 



a R (l/r) M ,N = -PV f X ' e- u/z {M/N} f {u)du 

Z J o 



j3 Q a R {q) 

477 



(24) 
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zeros 



poles 



[2/2] 
[3/3] 

[4/4] 
[5/5] 



[6/6] 



[7/7] 



[10/10] 



D 490^9 

U . t:Z UtJZ 


n 


. JUI I Ol 


n 

u . 


-UUUZtJ 






n 

u 


,JU4: / Ol 




Z : 


9^^4 




D Q^977 


n 

u 


74D9Qi 

. 1 4:UZc/l 


-0 


.78648 


fl 77QD8i 

U. I 1 JUOI 


D Q^977 


n 

u. 


74D9Qi 


-0 


.78648 




1 .UUOO 1 









.50170 




fl Q4714 






-1. 


.20630 




X . uuoou 


1 


1 D974i 

. 1 UZl I Til 


-0 


.90978 


1 DQ7Sni 
i.uy i oui 


-L . UUtJuU 


1 

-L . 


. lUi/ I Til 


-0 


.90978 


i.uy I oui 


9 SSD9Q 

Z . OOUZc/ 









.49992 




QQ840 






-2 


.50005 


U. OC/UZUl 


U . 1 tJUUO 


1 

J. . 


. OUOOOl 


-2. 


.50005 


U .OC/UZUl 


U. i ouuo 


1 

1 , 


.OUOOOl 


-1. 


.35858 


1 4fl47fii 

_L .4:Ut: 1 Ul 


1 97QQ^ 


n 


8401 Qi 

. Ot:U X C/l 


-1. 


.35858 


1 4fl47fii 

_L .4:Ut: 1 Ul 


1 97QQ^ 

-L . Z l C/C/O 


n 


R4D1Qi 

,Ot:U 1 C/l 


0.499997 




1 0009^ 






-1 


.72341 


U.OZiO t Ol 


n 8i n 

U.OIOIU 


1 

J. . 


. U I O I XI 


-1 


.72341 


U. OZtJ I Ol 


n hi n 

U.OltJIU 


1 

± . 


fi7^71 i 

.U 1 O 1 ±1 


-1 


.16083 


X .UZOUU1 








-1. 


.16083 


- 1 . UZOUU1 




n 

U. 


C/Z 1 OUI 


0. 


.50000 




1 ^1Q7 

1 .OO 1c/ ( 




-u. 


.c/z < oui 


1. 


.07610 




n QQQQzL 






-1. 


.77158 


u .yuuiui 


fl 74^1 7 


i 

i . 


c/OO i Ol 


-1 


.77158 


u. yuuiui 


74^1 7 


i 

± . 




-1 


.24094 


1 7Q4fi^i 

J. . 1 C/tiUOI 


1 ^1901 


i 

i . 




-1 


.24094 


1 7Q4fi^i 


1 ^1901 


i 

i . 


974.Q8; 





.50000 




l.oyl4z 


A 

u. 


ooozli 


1 


.60358 




1.59142 


-0 


,56621i 


76. 


.62431 




1.95927 






-2. 


.25821 




1.00000 






-2. 


.14754 


1.01218i 


0.91881 


2. 


53060i 


-2 


.14754 


-1.01218i 


0.91881 


-2 


53060i 


-1. 


.77230 


1.74984i 


1.54963 


1. 


,87556i 


-1. 


.77230 


-1.74984i 


1.54963 


-1. 


.87556i 


-1 


.16222 


2.4711H 


1.95905 


1. 


23582i 


-1 


.16222 


-2.47111i 


1.95905 


-1. 


.23582i 





.50000 




2.24957 


0. 


55669i 


1 


.50020 




2.24957 


-0. 


.55669i 


3 


.90255 





We see that PAs to f{u) actually reconstruct zeros and poles on the real axis 
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and that there are many other complex zeros and poles drifting away from u = 0, 
most probably "simulating" the essential singularity of f(u) at u — oo. 

Now we can see what values for (22) can one get from (|2Jj). We shall compare 
these values with the ones obtained in |L3| as the best partial sums (denoted by 
Q!fl,jv(l/ r ))) anc ^ those calculated from formula (0), though one must keep in mind, 
that (p~5f ) gives the exact value of otR in position space, while we are here estimating 
the sum ( p2j) which is fiLBj) " minus higher twists" . 

We present below the case M = N, but include also [2/3], as it is the lowest 
PA which gives already values quite close to the limit and uses only five coeffcients 
of the expansion. 



l/r[GeV] 


10 


20 


50 


100 


[2/2] 


.2903148980 


.2370475373 


.1814661174 


.1522256606 


[2/3] 


.2710750452 


.2232228306 


.1757241919 


.1497191873 


[3/3] 


.2714196149 


.2235681574 


.1759161959 


.1498234840 


[4/4] 


.2714564564 


.2234913019 


.1758530906 


.1497850862 


[5/5] 


.2714137805 


.2234817444 


.1758530534 


.1497859540 


[6/6] 


.2714097805 


.2234792391 


.1758723578 


.1497857339 


[7/7] 


.2713988440 


.2234777131 


.1758524051 


.1497857952 


[8/8] 


.2714007078 


.2234779029 


.1758523985 


.1497857906 


a R (l/r) 


.2714046 


.2218345 


.1746993 


.1491174 


OR,iv(l/ r ) 


.2856 


.22262 


.17647 


.14868 



l/r[GeV] 


10 3 


10 4 


10 5 


10 6 


[2/2] 


.09908701930 


.07409640462 


.05937036386 


.04957967229 


[2/3] 


.09908856884 


.07413913414 


.05938321829 


.04958332635 


[3/3] 


.09909570614 


.07413847053 


.05938281763 


.04958320088 


[4/4] 


.09909193228 


.07413842675 


.05938289532 


.04958322765 


[5/5] 


.09909214599 


.07413844576 


.05938289564 


.04958322732 


[6/6] 


.09909214737 


.07413844662 


.05938289573 


.04958322732 


[7/7] 


.09909214937 


.07413844659 


.05938289572 


.04958322732 


[8/8] 


.09909214928 


.07413844654 


.05938289571 


.04958322732 


«ij(V r ) 


.09901686 


.07413084 


.05938213 


.04958315 


«A,7v(l/ r ) 


.098937 


.07413536 


.05938222 


.04958339 



We see that olr{1/t)m,m converge when M grows, but not to the values of a# 
- we attribute this to the difference made by "higher twists". In other words we 



claim that the procedure we propose sums the series (]2"2" ) well, which allows us 



to see clearly the difference between the perturbative solution - represented by 
this series - and the full solution (^) including also nonperturbative effects. We 
demonstrate, therefore, once again the fact that the perturbation series contains 
only a part of the information about the full solution and that the requirement 
that the solution must be real for positive values of the coupling constant is not 
sufficient to compensate for the information contained in "higher twists" terms. 

4 Branch points in the Borel transform 

If the Borel transform has a branchpoint at z — K, instead of a pole, then the 
integral (|5|) for xe(0, oo) instead of (|7|) takes the form: 

/•OO /"OO 

$(x) = / e~ u Reg(ux)du±i e~ u Img(ux)du (25) 
Jo Jk/x 

if g{z)(z — K) is finite for z — > K (what is a condition for the existence of (||). 

Again, we see that the second integral, which is equal to the imaginary part of $(x), 
has all terms of the asymptotic expansion for x = 0+ vanishing, and therefore it is 
the first integral which is of interest for us. 

In the same way, as in the case of the simple pole, we can find an integral 
representation for the analytic function having the required asymptotic expansion 
(^) and being real on the positive semiaxis: 

roo+ie 1 roo 

£(z)= e- u g{uz)du- — e~ ulz \g{u + ie) - g{u - ie)}du (26) 

Jo+ie ZZ J K 

These facts were also observed earlier but we want to stress that ReQ(x) has 

the same asymptotic expansion at zero as $(x) and therefore it cannot be reliably 
estimated by calculation of partial sums. 

If we want now to use the same concept as before: to sum using PAs, we 
must remember that they will most probably not converge on semiaxis xe(K, oo). 
We conjecture, however, that if PAs to g(z) converge in measure arbitrarily close 
to positive semiaxis (where we expect a cut, "reconstructed" by PAs, should lie), 
then the integral: 

£(z) MN = PV e- u [M/NUuz)du (27) 
Jo 



will converge to ReQ(z) (see (P5|) ) for z^O. 

To check a plausibility of this hypothesis, we try to calculate this way: 

poo 

F(x) = Re / e- u/x log(l - u)du (2£ 



o 
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i.e. we approximate F(x) by: 

F m ,n{ x ) = Py r e- u l x u[M/N} lorAl _ t)/t {u)du (29) 

J 

As log(l — t)/t is the Stieltjes function, all poles of PAs to it lie on the semiaxis 
(l,oo) and interlace with zeros 01,0. 

On the other hand, F(x) can easily be converted to a simple principal value 
integral: 

poo p-u/x 

Fix) = -xPV / du (30) 

JO I — u 

We give below values of F M)N (x) for few values of x and compare them with 
values obtained from a numerical integration of the last formula. 



x 



.05 



.1 



.2 



.95 



[1/1] 
[1/2] 
[2/2] 

[2/3] 
[3/3] 
[3/4] 
[4/4] 
[5/5] 
[6/6] 
[7/7] 
[8/8] 
[9/9] 
[10/10] 

[11/H] 



0.0527903605 
0.0527972269 
0.0527977128 
0.0527977828 
0.0527977927 
0.0527977949 
0.0527977953 
0.0529797953 
0.0527977953 
0.0527977953 
0.0527977953 
0.0527977953 
0.0527977953 
0.0527977953 



0.1128659958 
0.1131155255 
0.1131547428 
0.1131540942 
0.1131486878 
0.1131464797 
0.1131466717 
0.1131470772 
0.1131470109 
0.1131470222 
0.1131470201 
0.1131470205 
0.1131470205 
0.1131470205 



0.2754741870 
0.2731042382 
0.2705616695 
0.2703267712 
0.2708042139 
0.2708502503 
0.2707469311 
0.2707743167 
0.2707637347 
0.2707667730 
0.2707662425 
0.2707662164 
0.2707662706 
0.2707662537 



0.4092156955 
0.8349862840 
0.7678390728 
0.6479952192 
0.7437635764 
0.7091822199 
0.6971934572 
0.7066673691 
0.6701782542 
0.7125737482 
0.7110753552 
0.7110508353 
0.7113684213 
0.7114624672 



0.049324951 
-1.544989018 
-1.080166535 
-0.590917162 
-1.199222773 
-0.256031716 
-0.961775112 
-0.711782921 
-0.586264408 
-0.577295644 
-0.621454057 
-0.665707635 
-0.688595754 
-0.691851149 



It sholud be noticed that starting from [8/8] the quadruple precision was nec- 
essary to obtain the above numbers. 



5 Conclusions 

We have proposed above to use PAs to sum the Borel transforms of the series 
appearing in QCD, and then to calculate Cauchy Principal Value of the Borel 
integral to obtain a function which has the given series as its asymptotic expansion 
at zero, and which is at the same time real on the positive semiaxis in the complex 



10 



coupling constant plane. The proposal is not entirely new M , but we have pointed 
out that it produces values of the well defined analytic function, one of the variety of 
them having the same asymtotic expansion. We have also argued that the method 
should work not only in the case when the Borel transform has poles on the positive 
semiaxis, but also when it has branchpoints there. This latter point is illustrated 
by the very simple numerical experiment. 

6 Acknowledgements 

The author is indebted to prof. M. Jezabek for drawing his attention to his paper 



13], and for interesting discussion and also to dr. P. Raczka for pointing him out 



paper fSfl. 

References 

[1] J. Zinn- Justin, Phys.Rep. 70 (1981) 109. 

[2] A.I.Veinshtein, V.I.Zakharov, Phys.Rev.Lett 73 (1994) 1207. 

[3] G.A. Baker Jr., P. Graves-Morris, "Pade Approximants" , Encyclopedia of 
Mathematics and its Applications vol. 13.14, Addison- Wesley Publ. Co. 1981 

[4] J. Gilewicz, "Approximants de Pade", Lecture Notes in Mathematics, Springer 
Verlag 1978 

[5] J. Nuttall, Journal of Math. Anal, and Appl. 31 (1970) 147 

[6] L. Lovitch, M.F. Marziani, II Nuovo Cimento 76A (1983) 615 

[7] M.F. Marziani, II Nuovo Cimento 99B (1987) 145 

[8] PA. Raczka, Phys. Rev. D 43 (1991) R9 

[9] C. Brezinski, Journal of Approx. Theory 54 (1988) 210 
[10] C.N. Lovett-Turner, C.J. Maxwell, Nucl. Phys. B432 (1994) 147 
[11] J. Ellis, E. Gardi, M. Karliner, M.A. Samuel, Phys. Lett. B366 (1996) 268 
[12] L.Richardson, Phys. Lett. B82 (1979) 272 

[13] M. Jezabek, M.Peter, Y. Sumino, hep-ph /9803337. March 1998, Phys.Lett. 
B428 (1998) 352 

[14] G. Grunberg, Phys. Lett. B325 (1994) 441 



11 



